Effect of the Heterogeneity of Metamaterials on Casimir-Lifshitz Interaction 
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The Casimir-Lifshitz interaction between metamaterials is studied using a model that takes into 
account the structural heterogeneity of the dielectric and magnetic properties of the bodies. A re- 
cently developed perturbation theory for the Casimir-Lifshitz interaction between arbitrary material 
bodies is generalized to include non-uniform magnetic permeability profiles, and used to study the 
interaction between the magneto-dielectric heterostructures within the leading order. The metama- 
terials are modeled as two dimensional arrays of domains with varying permittivity and permeability. 
In the case of two semi-infinite bodies with fiat boundaries, the patterned structure of the mate- 
rial properties is found to cause the normal Casimir-Lifshitz force to develop an oscillatory behavior 
when the distance between the two bodies is comparable to the wavelength of the patterned features 
in the metamaterials. The non-uniformity also leads to the emergence of lateral Casimir-Lifshitz 
forces, which tend to strengthen as the gap size becomes smaller. Our results suggest that the recent 
studies on Casimir-Lifshitz forces between metamaterials, which have been performed with the aim 
of examining the possibility of observing the repulsive force, should be revisited to include the effect 
of the patterned structure at the wavelength of several hundred nanometers that coincides with the 
relevant gap size in the experiments. 

PACS numbers: 05.40.-a, 81.07.-b, 03.70.+k, 77.22.-d 



I. INTRODUCTION 



Despite nearly six decades of research after the original 
works of Casimir |l| and Lifshitz [2|, the dependence of 
Casimir-Lifshitz force between bodies on their geometri- 
cal and material properties is still a subject of ongoing 
investigation |3[ . The effect of geometry has been studied 
using a variety of techniques, which include perturbative 
expansion around ideal geometries (3411 and in dielectric 
contrast @-Qj]], semiclassical 12 1 and classical ray-optics 
[lH approximations, multiple scattering [3 El I and mul- 
tipole expansions |16l - [l8| , world- line method 19j . exact 
numerical diagonalization methods (2fjj , and the method 
of numerical calculation of the Green function [2l[ . These 
studies have significantly advanced our understanding of 
the subtle effect of geometry on Casimir-Lifshitz interac- 
tions, and have led to proposals for using the knowledge 
in designing useful nano-scale mechanical devices [22j . 



The dependence on material properties has also been 
studied extensively since the work of Dzyaloshinskii, Lif- 
shitz, and Pitaevskii, who pointed out that the force can 
be attractive or repulsive depending on the relative val- 
ues of the dielectric constants of the successive layers 
[23| . The existence of a repulsive mode of the interac- 
tion is very interesting, as it explains, for example, why 
a wetting layer of liquid should form on a solid in equilib- 
rium with vapor [23j . Despite the theoretical possibility, 
it is not trivial to find a condition where the Casimir- 
Lifshitz interaction between two solid bodies that are ei- 
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ther metallic or dielectric turns repulsive due to the pres- 
ence of a non-solid medium between them. However, a 
recent experiment has shown that such a repulsive force 
can be observed between gold and silica particles that are 
separated by bromobenzene [24] . Another possibility for 
repulsive Casimir-Lifshitz interactions was pointed out 
by Boyer, who showed that the force between a purely 
dielectric semi-infinite body and a purely magnetic one 
separated by vacuum, is repulsive [25| ■ Considering that 
the Casimir-Lifshitz force will be a key player in the 
realm of micro / nano-electromechanical systems (the do- 
main that involves length scales of the order of 100 nm to 
1 /im) the possibility of producing repulsive forces gives 
hope for eliminating stiction. However, from the work 
of Lifshitz we know that at those distances the Casimir- 
Lifshitz force will be determined by the relatively high 
frequency part of the permittivity and permeability of 
the materials in imaginary frequency, and at such high 
frequencies the response of natural magnetic materials to 
the electromagnetic field is neglig ible (n ~ 1) HlUm. 
In other words, while the theoretical possibility for cre- 
ating a repulsive force exists, natural materials with the 
required magnetic properties cannot be found. 

In recent years, engineered materials — called 
metamaterials — have been developed based on the 
proposed concept of negative refractive index (27[, and 
their physical properties have been extensively studied 
|28l - l31| . This development has brought about the pos- 
sibility of designing materials with special permittivity 
and permeability in a desirable range of frequencies. 
This could, in turn, result in producing nontrivial 
magnetic response in a broad range of frequencies, and 
possibly help achieve the repulsive Casimir-Lifshitz 
force |32h37| . A main characteristic of metamaterials is 
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FIG. 1: Two semi-infinite metamaterials modeled as 
chessboard-patterned magneto-dielectric media. The struc- 
ture can be characterized with the separation H, wavelengths 
\ x and X y , fractions f x and f y , permittivities ei (uj) and 
€2 (ui), and permeabilities fj,± (uj) and fi2 (w). The vector 
(aX x ,bX y ) denotes the displacement of the upper object rela- 
tive to the lower one. 



their periodic engineered structure, which could involve 
features at length scales ranging from hundreds of 
nanometers to a few microns [31(. These features could 
correspond to metallic split-ring resonators that are im- 
planted in a dielectric background, or similar structures 
present in photonic crystals. This means that while the 
macroscopic response of metamaterials to electromag- 
netic fields can be described via appropriate frequ ency 
dependent permittivity and permeability functions [28], 
at shorter distances they should be treated as a periodic 
distribution of regions with contrasted permittivity and 
permeability response functions. This is of particular 
importance in the calculation of Casimir-Lifshitz force, 
as we know that any lateral feature in the material 
properties will affect the force when the gap size is of 
the order of the characteristic length scale set by the 
heterogeneity; as argued above the typical length scales 
for these features coincide with the range at which 
Casimir-Lifshitz forces are most significant. 

The Casimir-Lifshitz force between metamaterials has 
been investigated recently [321 - 1371 ] . In these studies, the 
material properties are taken into consideration at the 
macroscopic level, in the sense that the permittivity and 
permeability corresponding to uniform materials have 
been incorporated in Lifshitz theory. In this paper, we 
examine the effect of the periodic structure of metama- 
terials on the Casimir-Lifshitz interaction using a gen- 
eralization of the dielectric contrast perturbation theory 




FIG. 2: Schematics of the assortment of magneto-dielectric 
objects described by their different dielectric function and 
magnetic permeability profiles. 



8|. |38l| and its application to dielectric heterostructures 
33| . We develop a perturbative scheme for the calcula- 
tion of the Casimir-Lifshitz force as a series expansion in 
powers of the contrast in permittivity and permeability 
profiles. We use the theoretical formulation to calculate 
the force for a model of metamaterials that is made of a 
two dimensional periodic structure of varying magneto- 
dielectric properties, as shown in Fig. [TJ We find that 
the periodicity in the structure of our model metamateri- 
als causes the normal Casimir-Lifshitz force to change as 
compared to its value when the materials are assumed to 
be uniform. This change is found to be significant when 
the distance between the two bodies is comparable to 
the wavelength of periodic structure of the bodies. The 
heterogeneity also introduces a lateral component to the 
Casimir-Lifshitz force, which is analogous to the lateral 
Casimir force between co rrug ated surfaces [U, H(| and di- 
electric heterostructures }39j . and is more significant at 
smaller separations. 

The rest of the paper is organized as follows. In 
Sec. HU we develop the theoretical formulation of 
the perturbation theory that can be used for studying 
Casimir-Lifshitz interaction between magneto-dielectric 
heterostructures. Section IIIII is devoted to applying the 
perturbative scheme to the particular problem of semi- 
infinite patterned magneto-dielectric structures at the 
leading order of the perturbation theory. In Sec. IIV1 
the results of the calculation of the normal and lateral 
Casimir-Lifshitz forces are shown, and finally, Sec. IVl 
concludes the paper with some discussion and remarks. 



II. THEORETICAL FORMULATION 

We consider an arrangement of magneto-dielectric ob- 
jects in space with arbitrary shapes and frequency depen- 
dent dielectric and magnetic properties, as shown in Fig. 
[2] The frequency- and space-dependent dielectric func- 
tion e(w, r) and magnetic permeability r) describe 
the medium. As shown in Ref. [38[ , the Casimir-Lifshitz 
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energy of the system can be written as 
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involves the dielectric function and magnetic permeabil- 
ity profiles in imaginary frequency. 



A. Perturbation Theory 

Following [sl, HH], we develop a systematic expansion 
of Eql in terms of the dielectric contrast Se(i(, r) = 
e(«C,r) — 1 and the inverse magnetic permeability con- 
trast («C,r) = (i(,r) — 1. Using the Fourier 
transforms 

MC;q,q') = / d 3 rd 3 r'/C y (C;r,rV q V q ' r \ (3) 
Se(i(,q) = J d?v [e«,r)-l] e^ r , (4) 
*A _1 (*C,q) = / d 3 r [^(iC.rJ-l] e^ r , (5) 
the kernel in Eq. © can be decomposed as 

£ij(C;q,q') = /Co, l3 (C,q)(2^) 3 5 3 (q + q') 



5/C?.(C;q,q')+^™(C;q,q')- (6) 



Here 



c 2 

&0,ij (C) q) = -o <*»j + 9 2(5 y ^ 9i«7 > ( 7 ) 



corresponds to the empty space, and 

<5/c^(C;q, q ') - ^M?«,q + q'), (8) 

J/C™(C;q,q') = (^-gfc^yJtf/i-^iCq+qO.CS) 

entail the permittivity and permeability profile. 

We can now recast the expression of Eql into a per- 
turbative series using the identity 
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The general term for the series expansion in Eq. ([TO 
takes on the form 
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Using the above explicit form, the Casimir-Lifshitz energy can be calculated for an arbitrary assortment of magneto- 
dielectric materials by following standard diagrammatic methods. 

In the rest of this paper, we focus only on the second order and calculate the energy for periodic structures. 



B. Second Order Term 



We consider the leading contribution in the perturbation theory, which comes from the second-order term of the 
series (or its correction by a so-called Clausius-Mossotti factor using a resummation [H, HH). We note that at the 
second order, we can use Sjl^ 1 — —5/1, which will simplify the calculations that will follow later on. We find the 
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second order Casimir-Lifshitz energy as 
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The above result strongly depends on the relative posi- 
tioning, geometry, and electromagnetic characteristics of 
the bodies. Note that the Sidp, term manifests the non- 
additive nature of electric and magnetic contributions to 
the Casimir-Lifshitz energy. 



III. MAGNETO-DIELECTRIC 
HETEROSTRUCTURES 

We now consider two parallel semi-infinite magneto- 
dielectric bodies placed at a separation H, such as the 
one shown in Fig. [TJ Introducing the labels u and d for 
"up" and "down" bodies, the permittivity function e and 
the permeability function /i can be written as 



e(i(,r) 



/i«,r) = < 
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where r = (x, z). 
Using the identity 
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and introducing the new variables P = i(k — q) and 
Q = k + q, we can simplify Eq. (1131) by performing the 
Gaussian integrations over the variable P. This yields 
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This result can now be used to study the Casimir-Lifshitz 
interaction between two macroscopic bodies with any 
magneto-dielectric profiles at the second order in per- 
turbation theory. 



A. Two homogenous semi-infinite bodies 

Let us first consider two homogenous semi-infinite bod- 
ies of area A and separation H. Assuming that the per- 
mittivity and permeability are frequency- independent, 
the Casimir-Lifshitz force F — —dE-ijdFl takes on a sim- 
ple form 



F 



hcA 



6407r 2 iJ 4 
x [23(6e d 5e u - 



(20) 



as found previously in the literature |4jj-|43| . We can now 
consider patterned magneto-dielectric objects, and study 
how the structural heterogeneity affects both the normal 
and lateral Casimir forces. 



B. Patterned magneto-dielectric objects 

Let us now consider the configuration shown in Fig. [TJ 
The magneto-dielectric "chessboard" heterostructure can 
be characterized with wavelengths X x and X y along the x 
and y directions. In a repeat unit of the material along 
each direction a (a = x,y), a fraction f a of the material 
has permittivity £2 (w) and permeability fi2 (w), and the 
remaining fraction (1 — f a ) has permittivity e\ (lu) and 
permeability /zi (w). Below, we will consider two possi- 
bilities; one with the domains with higher permittivity 
and permeability coinciding and another where they are 
in a staggered configuration. The vector (aX x ,bX y ) de- 
notes the displacement of the upper object relative to the 
lower one, as shown in Fig. [TJ 

We use the Clausius-Mossotti resummation of the per- 
turbation theory for the permittivity contribution [8L [38lj . 
which amounts to replacing 5e(i£, x) in Eq. (fT5|) by 
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FIG. 3: (Color online). The normal Casimir force (in units 
of F c ) for the case of EhMh-ElMl, as a function of the 
dimensionless parameter a for various distances H. Here 
Xx — \ — 500 nm, and the two panels correspond to (a) 
fx = f y = 0.5, and (b) f x = 0.75 and /„ = 0.25. 



FIG. 4: (Color online). The normal Casimir force (in units 
of Fc) for the case of ElMh-EhML as a function of the 
dimensionless parameter a for various distances H. Here 
\x — Xy — 500 nm, and the two panels correspond to (a) 
fx=fy = 0.5, and (b) f x = 0.75 and f y = 0.25. 



Due to the periodicity of the magneto-dielectric struc- 
tures, it is natural to use the Fourier series expansion of 
the permittivity and permeability profiles. We have 
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The Fourier series coefficients can be easily found as 
AimK) = 2 |_0e 2 «) - dei(zC)J — 



and 
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(22) 

for n / and to 7^ 0. Using Eq. (fTo| . one obtains the 
Casimir-Lifshitz energy .Ec of the chessboard magneto- 



dielectric heterostructure as 
HA 



E c = 
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2(2p 2 - l)A nm (i() S„ m «)] , (23) 



where the prime on the summation indicates that the 
to = n = term comes with a pref actor of 1/2. 



C. Material Properties 

While the experimental realizations of metamateri- 
als involve a multitude of complex structures, in the 
present study we consider a simplified model where 
the heterostructure could have two types of effec- 
tive permittivities — corresponding to metallic and di- 
electric materials — and two types of permeabilities — 
corresponding to magnetic and non-magnetic materials. 
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FIG. 5: (Color online). The lateral Casimir force for the case 
of EhMh-ElMl, as a function of the dimensionless parameter 
a for various distances H. Here = A y = 500 nm, and the 
two panels correspond to (a) f x — f y = 0.5, and (b) f x — 0.75 
and /„ = 0.25. 



plat 

— (Nm- 2 ) 







H (nm) 






— 100 






120 






150 


\ \ 




200 



























-0.4 -0.2 0.0 0.2 0.4 

(b) a 

FIG. 6: (Color online). The lateral Casimir force for the case 
of ElMh-EhMl, as a function of the dimensionless parameter 
a for various distances H. Here \ x = \ y = 500 nm, and the 
two panels correspond to (a) f x = f y — 0.5, and (b) f x = 0.75 
and /„ = 0.25. 



For the metal, which has a relatively higher dielectric con- 
stant especially at lower frequencies, we use the Drude 
model that in imaginary frequency reads 



n 2 

e„(i() = l+ " 



(24) 



Alternatively, for the dielectric medium with relatively 
lower dielectric function we use the Drude-Lorentz model 



ez(*<) = 1 



fX 2 



(25) 



Similarly, we choose a simple Drude-Lorentz model for 
the permeability of the magnetic material, namely 



Mfc(*C) = l + 



(26) 



C 2 + < + 7m c' 

while for the non-magnetic material we have 

W(*'C) = 1- (27) 

In the above equations, w e (ui m ) is the electric 
(magnetic) resonance frequency, and 7 e (7 m ) is the 



electric (magnetic) dissipation parameter. Using the 
plasma frequency of gold u; p (Au) = u p = 1.37 x 
10 16 rad/s as a frequency scale, the numerical val- 
ues of the magneto-dielectric characteristic param- 
eters are chosen as: Q,d/ui p = 1.0, Jd/u p = 0.004, 
Q e /ujp — 0.04, n m /oj p — 0.1, <jJ e /u p = uj m /(jj p = 0.1, and 
j e /uj p = j m /uj p = 0.005 [13. 



We consider two different possibilities: (1) When the 
metallic patch has magnetic properties and the dielectric 
patch is non-magnetic. In this case, which we represent it 
schematically as EhMh-ElMl, we have £2 = eh, A* 2 = 
ei = ej, and /zi = (J,i- (2) When the dielectric patch 
has magnetic properties and the metallic patch is non- 
magnetic. In this case, which we represent it schemati- 
cally as ElMh-EhMl, we have 62 = ei, fJ-2 = fJ-h, £1 = 
and Hi = hi. Below, we will study both the normal and 
lateral Casimir-Lifshitz forces in both of these cases. 
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FIG. 7: (Color online). Vector field plots of the lateral Casimir-Lifshitz force as a function of the relative displacements a and 
b, along the x and y directions respectively, for the EhMh-ElMl model at H — 100 nm. (a) f x — fy — 0.5 and X x = X y = 500 
ran. (b) f x = 0.75 and /„ = 0.25 and \ x = X v = 500 nm. (c) f x = 0.75, /„ = 0.25, \ x = 500 nm, and \ v = 200 nm. 



IV. CASIMIR-LIFSHITZ FORCES BETWEEN 
THE CHESSBOARD STRUCTURES 

Using the material properties described in Sec. MI CI 
above, we can now calculate the Casimir-Lifshitz en- 
ergy for the chessboard magneto-dielectric heterostruc- 
ture shown in Fig. [T]for the two cases denoted as EhMh- 
ElMl and ElMh-EhMl above. Due to the lateral hetero- 
geneity of the magneto-dielectric properties, the two bod- 
ies exert both normal and lateral Casimir-Lifshitz forces 
on each other. 



A. Normal Force 

The normal force between the two patterned structures 
can be calculated as Fq™ — — ^§-. We set 6 = 0, and 
study the normal force as a function of a for various val- 
ues of the gap size H . We choose to normalize the force 
using Fq, which is the contribution of the m = n = 
term in Eq. (|23|) to the normal force. Figure [3] shows the 
normal force relative to F®, for the EhMh-ElMl case. Fig- 
ure^ corresponds to the symmetric configuration where 
fx = fy = 0.5 whereas Fig. |3b corresponds to an asym- 
metric configuration where f x = 0.75 and f y = 0.25, 
and in both cases X x — X y = 500 nm. Figure [3] shows 
that depending on how the different patches with dif- 
ferent magneto-dielectric properties are positioned with 
respect to one another, the normal Casimir-Lifshitz force 
can change in magnitude relative to the case with uni- 
form magneto-dielectric configuration. When patches 
with similar properties are opposite one another (a = 0) 
the attractive normal force is at its maximum, while the 
force is weakest when dissimilar patches are exactly op- 
posite one another. While this relative change depends 
strongly on the gap size, we note that it can easily amount 
to a few percent in the experimentally relevant gaps sizes 
of a few hundred nanometers, and could even reach the 



value of 35% for the gap size of H = 100 nm for the 
asymmetric example. Figure U shows a similar behavior 
for the ElMh-EhMl case, which shows a relatively less 
dramatic change in the asymmetric example. 

B. Lateral Force 

The lateral force between the two patterned structures 
is a vector, with its value and direction depending on the 
relative positioning of the two bodies. For simplicity, we 
first focus on the case with 6 = 0, and only study the 
lateral force for unidirectional displacements along the x 
axis (see Fig. [T]). In this case, the force also lies along the 
x axis for symmetry reasons, and we can find its value 
using F ut = — -r-^ 2 -. as a function of a for different 

G A x da J 

gap sizes H. Figure [5] shows the lateral force per unit 
area in SI units, for the EhMh-ElMl case. Figure 
corresponds to the symmetric configuration where f x — 
f y = 0.5 whereas Fig. [SJd corresponds to an asymmetric 
configuration where f x — 0.75 and f y = 0.25, and in 
both cases \ x = X y = 500 nm. Figure [5] shows that the 
lateral Casimir-Lifshitz force is very sensitive to the value 
of and its dependence on the lateral displacement 
reflects the symmetry or asymmetry of the relative sizes 
of the two patches. While the form of the lateral force at 
relatively larger gap sizes tends to a sinusoidal form, at 
smaller separations higher harmonics contribute as well 
to reflect more of the details of the heterogeneity. Figure 
|6] shows the lateral force for the ElMh-EhMl case, which 
shows a similar behavior as compared to the previous 
case. 

In Fig. [71 the vector field for the lateral Casimir- 
Lifshitz force, defined as F lat (a,6) = -^^x - 

j-^§^y, is plotted as a function of a and 6. One can 
see that the symmetry of the heterostructure affects the 
configuration of the lateral Casimir-Lifshitz force as a 
vector field, and that the numerous parameters involved 



8 



can provide opportunities for a rich variety of engineered 
patterns for the lateral force. 

V. DISCUSSION 

We have studied the Casimir-Lifshitz interaction be- 
tween two metamaterials modeled as periodic arrayed 
structures containing domains of varying magneto- 
dielectric properties. We have considered two types 
of permittivity functions — metallic and dielectric — and 
two types of permeability functions — magnetic and non- 
magnetic — and their corresponding two combinations. 
For both combinations, we have found significant changes 
in the value of the normal Casimir-Lifshitz force relative 
to the value that corresponds to the uniform (macro- 
scopic) model of the materials. The relative change is 
increased as the gaps size is decreased, and reaches a few 
percent in the experimentally relevant gaps sizes of a few 
hundred nanometers, while it could even reach 35% when 
the gap size is 100 nm in the models we studied. Consid- 
ering how delicate it is to find the condition to achieve 
the repulsive force in realistic situations as recent studies 
have revealed [lij , our results show that the effect of the 
structural heterogeneity should be taken into account in 
determining whether the force is repulsive or not. This is 
particularly pertinent as the characteristic length scale of 
the periodic features in realistic metamaterials coincides 
with the gap sizes at which the Casimir-Lifshitz force is 
particularly relevant. While this issue has been ignored 



by all previous studies, we note that our study has been 
performed within the scope of the magneto-dielectric con- 
trast perturbation theory at its leading order and should 
be considered more as an indication of the relative signifi- 
cance of the effect rather than a study that could provide 
numerically accurate results [44( . To that end, one needs 
to employ more sophisticated numerical methods simi- 
lar to those developed to study the effect of geometry 

dam. 

Another consequence of the structural heterogeneity 
of the metamaterials is the possibility of the emergence 
of lateral Casimir-Lifshitz forces, which are stronger for 
smaller gap sizes. While all previous works on lateral 
Casimir force have focused only on unidirectional geo- 
metrical or material heterostructure features, we have 
considered a two dimensional pattern and presented the 
vector field distribution of the lateral force. These forces 
are very sensitive to the details of the periodic patterns 
of the magneto-dielectric properties, and their versatility 
allows them to be amenable to detailed engineering by 
changing these features. 
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